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We investigate the primordial power spectrum of the density perturbations based on the assump- 
tion that spacetime is noncommutative in the early stage of inflation. Due to the spacetime non- 
commutativity, the primordial power spectrum can lose rotational invariance. Using the k-inflation 
model and slow-roll approximation, we show that the deviation from rotational invariance of the 
primordial power spectrum depends on the size of noncommutative length scale Ls but not on sound 
speed. We constrain the contributions from the spacetime noncommutativity to the covariance ma- 
trix for the harmonic coefficients of the CMB anisotropics using five-year WMAP CMB maps. We 
find that the upper bound for Ls depends on the product of sound speed and slow-roll parameter. 
Estimating this product using cosmological parameters from the five-year WMAP results, the upper 
bound for Ls is estimated to be less than 10~^^cm at 99.7% confidence level. 



in 
O 

u 

^ ■ 

Oh 
I 

O ' 

(N ■ 
> • 

^ : 

o . 

oo • 

(N : 

rn 
O , 

a\ . 
o ■ 



I. INTRODUCTION 

The inflationary cosmology [l| is the scenario of the 
very early universe. It provides a successful mechanism 
for generating nearly scale invariant primordial density 
perturbations, that give rise to galaxy formation and 
temperature anisotropics in the CMB which are in agree- 
ment with observation ■ If the period of inflation is suf- 
ficiently longer than that required for solving the horizon 
and flatness problems, such that the wavelengths of per- 
turbations which are observed today emerged from the 
Planck regime in the early stages of inflation, the physics 
on trans-Planckian scales should leave an imprint on the 
primordial density perturbations 

Near the Planck scale, the properties of spacetime are 
expected to be modified due to the quantum nature of 
gravity [3]. It has been shown that a consequence of 
string theory which is a promising candidate of quantum 
gravity, is that the spacetime is noncommutative Q 



Gaussianity of the CMB fiuctuations is assumed. Some of 
the estimators for quantifying the statistical anisotropy 
contributions in the CMB fluctuations have been pro- 
posed in [l^, [13, HBl • According to [13, [ll| , the statistics 
of the observed CMB fluctuations does not deviate from 
isotropy signiflcantly. 

In this work, we constrain the contributions from 
spacetime noncommutativity to CMB temperature fluc- 
tuations using five-year WMAP CMB maps. 



II. THE CONTRIBUTIONS FROM SPACETIME 
NONCOMMUTATIVITY 

In this section, we investigate the contributions from 
spacetime noncommutativity to the primordial power 
spectrum in the k-inflation model [l6j. 

We start with the general action of the inflaton of the 
form 



(1) 



2F(X, </))], 



(2) 



X 
S3 



where is an antisymmetric tensor. 

The influences of spacetime noncommutativity on the 
feature of power spectrum of primordial fluctuations have 
been studied by many authors [1, [3l- If but 

the primordial power spectrum can be- 
come direction-dependent, and consequently the statis- 
tics of CMB fluctuations becomes anisotropic. We are 
interested in this spacetime noncommutativity induced 
statistical anisotropy. 

Usually, the statistics of the CMB temperature fluctua- 
tions is supposed to be isotropic. However, recently there 
are many attempts to check whether the statistics of the 
CMB fluctuations is perfectly isotropic by searching for 
the statistical anisotropy contributions in the CMB sky 
maps [13, [m [13] ■ In the case where the statistics of the 
CMB fluctuations is anisotropic, the angular power spec- 
trum does not contain all the information about the sta- 
tistical properties of the CMB fluctuations even when the 



where (j) is the inflaton field and X — — (1/2)9^09^ 
Here, we have set the reduced Planck mass (SttG)"^/^ — 
1. To study the evolution of density perturbations during 
inflation, one expands the action ([2]) around the homoge- 
neous and isotropic background. In our consideration, we 
use the ADM metric formalism in which the line element 
is given by [i, [13] 



+ N'dt){dx^ + N^dt). (3) 



In the ADM formulation hij and (j) are the dynamical 
variables, while N and N'^ are Lagrange multipliers. To 
compute the perturbed action in the slow-roll approx- 
imation, it is convenient to use the uniform curvature 
gauge in which [TtI 5(j) = f{t,x) and hij = a^Sij, where 
a is the cosmic scale factor, S<j){t,x) = (j){t,x.) — 0o(i) is 
the perturbation in the inflaton field and the subscript 
represents the background value. 
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In order to find the second order action for (/s, one first 
computes the constraint equations for N and N"^ from the 
action ([2]) and solve these equations for N and to first 
order of (/3 Substituting the rcsuhs back in the action, 
expanding the action to the second order of perturbation 
and keeping the fowest order of sfow-roU parameter, the 
second order action for perturbation is given by 

/3 
d'x^ [PxoXoi^ovf - Pxod.^d^^) ,(4) 

= J d'a:p.xA^ ~ c^d^r) , (5) 

where = P,x / P.x is the sound speed, the subscript , X 
denotes a derivative with respect to X, p is the energy 
density of the infiaton and p^x = P,x + '^P,xxX. In the 
above calculation, we neglect the terms in the perturbed 
action that are multiplied by the metric perturbation be- 
cause these terms are subleading in slow-roll parameter 

We now study how the spacetime noncummutativity 
influences the action for the field perturbations. In order 
to take the effect of spacetime noncommutativity into 
account, we replace the ordinary products in the action 
with the star products. In curved spacetime, the star 
product can be expanded as Q 

{D^,---D^J){D,,---D,,g). (6) 

where Z?^ is the covariant derivative. In our considera- 
tion, we suppose that the nonzero components of Q'^" are 
= — O^^ = Ll/a?^ where Ls is the noncommutative 
length scale Q. 

Since the terms in the second order action that are mul- 
tiplied by metric perturbation are subleading in slow-roll 
parameter, the spacetime noncommutativity effect can 
be incorporated by replacing the products between fields 
in the Lagrangian by the star products, and expanding 
the Lagrangian to the second order in the field perturba- 
tion. In the calculation of star product, we also ignore 
the metric perturbation because it gives rise to the terms 
that are subleading in slow-roll parameter. Since it is 
possible to study the effect of spacetime noncommuta- 
tivity perturbatively [§1 , we do the calculation up to the 
lowest non trivial order in 8'"^. Up to the second order 
in 0'^", the star products of the multiple functions can 
be written as Q 

^ a<b 

^ D-DtD';,DAh---f„, (7) 

a<b,c<d ' 

where a, 6, c, d run over 1, . . . , n. 

For illustration, we suppose that P(X, 0) can be writ- 
ten as P(X,0) = F{(\))G{X) = ^"X'", and use eq. 



to compute the star product between field varibles. We 
first consider the case where the derivatives in eq. ([7]) act 
on the terms in G(X) only. In this case, the contribu- 
tion from spacetime noncommutativity starts to appear 
at the second order of 6'"', i.e., 

a<b,c<id 

= -^e^-^e''^(^ "'^'^~^h ',--'D,D,x,D^D,x, 

+m{m - l)X^-^Xi{-D^Dpd^^D,D„d^cl>„)^ , 
= -^e'"'&'"'(^lp.x„XoD^,DpXiD,D^Xi 

+PxoX„Xi{^Df,Dpd^^D,D„d^(bo)) , (8) 



where Xi = —do4>od''^ip is the firse order perturbation in 
X. In this equation, we omit the terms that are pro- 
portional to DpDpd-y(t>QDi,Dcrd^ 4)[) because these terms 
vanish after integration by parts. Using similar consider- 
ation, one can show that if the derivatives in eq. ([7]) act 
on the terms in Fi^) only, the noncommutative contribu- 
tion will be proportional to P^q^q, and if the derivatives 
act on the terms in both F{(f)) and G{X)^ the noncom- 
mutative contribution will be proportional to P^qXo- In 
our case, both contributions can be neglected because 
the terms that are proportional to P^o^o and P^qJCq are 
subleading in slow-roll parameter. Hence, the spacetime 
noncommutativity modifies the action in leading order of 
slow-roll parameter as 

where SS^'^'^ is the action for the perturbed field in eq. ([5]) 
and 5qS^^ = J d^xa^SeG. The action SqS^^ can be spit 

(2) 

into three parts as SqSx = SsSxi + SqSx2 + SsSxs, 
where the subscripts 1, 2, 3 denote the contribution from 
the first, second and third terms on the RHS of eq. ([8]). 
The expression for SqG in eq. ([5]) is valid for a generic 
P{X,(j)), because one always do Taylor's expansion of 
P{X, (j)) around the background when applying the star 
product in the action for perturbation. Moreover, the 
action ^ can be obtained by just replacing products 
with star products in the action 

We suppose that DpQ^'^ = 0, and do integration by 

parts for the action SqS^-^'' in the similar way as After 
doing integration by parts and evalulating the covariant 
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derivatives, each part of this action becomes 



the simple form 



SeSx3 



8 

aH^Li 



16 



-2Hdpipdpip — H^dpipdpip 



(10) 



(11) 



d^x^^^Pxx2Xdpdoip{dpd'>^ + Hdp^), 

o 



where p — 1,2 and H = a/a\s the Hubble parameter. In 
the above three equations and following calculation, we 
omit the subscript for P, p and X because from now 
on we will treat them as the background quantities. We 
note that in the calculation of SqSxs we suppose that 
X < X due to the slow-roll approximation. From these 
results, we get 



<5e4'^ = ^ / d^xaH^Lip^x 



(3dpd^vdpd\ (12) 
cldpdiipdpd^ip + 2l3Hdpipdpip — clH^dpipdpip 



where /? = 2c^ — 1. In order to find the evolution equa- 
tion and the corresponding solution from action ([9]), it 
is convenient to write the action in terms of the vari- 
able V — z(, where z ~ a^j2Xpx/ H [16]. The vari- 
able C is the curvature perturbation in comoving gauge 
which has a direct connection with the generation of large 
scale structure and CMB fluctuations. This quantity is 
related to ^ through ^ = Hip/(j)o. Expressing the ac- 
tion © in terms of v and expanding v in Fourier space 
as v{t,x) = /(d3k/(27r)3/2) (akWk(Oe''' '' + ^-C.), where 
at, aj^ satisfy the canonical commutation relations, the 
action ([9|), evaluated in the vacuum after normal order- 
ing, becomes 



drjd^k 
8a2 



.' |2 



(3K\'~clk'\vu\' 



(13) 



+ ( /3 ( 6— - 10{Hay 



/3-1 



(Hay \vk 



where the prime denotes derivative with respect to the 
conformal time rj = J da /a and k']_ = k\ -\- k"^. In 
the above action, we have neglected the terms c'^^ and 
p^x because these terms are subleading in slow-roll pa- 
rameter. Up to the lowest order of slow-roll parame- 
ter, we can write a? — \/{Hri)'^ + ©(e) — l/(i?ry)'^ and 
a" /a ~ 2/r], where e — —H/H^ is the slow-roll parame- 
ter. Using these approximations and writing in terms 
of the new variable j/k(?y) = {1 — PK^rf'k\Y/'^v\f^{vi), where 
= if^L^/S, the action up to first order in takes 



S = / dTyd^k 



(14) 



In the derivation of this action, we suppose that 
K}k\ff' < 1 and |c^| do not differ much from unity, so 
that one can expand 1/(1— /^K^fc^yy^) ~ l+PK'^kj_ri'^ + . . . 
and neglect the terms that are proportional to Pn^k^rf' . 
This implies that this action is valid only if < rjf = 
8/{H^L*k'^). For a suitable choice of initial conditions, 
we have — oo < i] < during inflation, so that the per- 
turbation mode k exits the horizon and sound horizon 
at rjc = —1/k and rjcs = — l/(fccs) respectively. Hence, 
we will be able to follow the evolution of the perturba- 
tion mode k well before its horizon (sound horizon) exit, 
i.e., rii < rjc {rji < rjcs), if the noncommutative length 
scale is required to be smaller than the Hubble radius, 
Ls < H-\ 

In order to study the influence of spacetime noncom- 
mutativity on the behavior of density perturbation, we 
compute the power spectrum of the curvature perturba- 
tion C from action From this action, one can show 
that the evolution equation for is 

y'i + [clkW '■^)v^^ 0' (15) 

where ^"^ = \ + k'^siv?{9) and 6 — sin^^(A:/A:^) denotes 
the angle between the vectors k and kj,. Following stan- 
dard procedure, the solution of eq. (fT^l) is given by ^185 



1 



y/2kc7f V kcsjrj 
so that the primordial power spectrum is givven by 



(16) 



vt 



2tt' 



k^ \vk/z\" 
27r2 



1 



(17) 

where — G is the Planck mass. It can be seen 
that the obtained power spectrum is direction-dependent 
due to the spacetime noncommutativity. The magni- 
tude of the noncommutative contribution depends on 
the angle between the vectors k and fcs, but not on 
the sound speed. From the expression for k, it can 
be seen that the deviation from rotational invariance of 
the power spectrum also depends on the noncommuta- 
tive length scale and Hubble parameter during inflation. 
The amplitude of this power spectrum can be defined as 
As — / {TTmpCse) . This quantity as well as k are eval- 
uated when the perturbation mode k exits the horizon. 
If we neglect the contributions from spacetime noncom- 
mutativity, the above power spectrum gives rise to the 
usual power spectrum as in [l6|. 
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III. THE CMB CONSTRAINTS 

In this section, we will constrain the contributions from 
spacetime noncommutativity to CMB anisotropics using 
the CMB data. It is well known that if the primordial 
power spectrum is direction-dependent, the statistics of 
CMB will become anisotropic, i.e., the two-point func- 
tion of the temperature fluctuations is no longer rota- 
tionally invariant. Here, we consider only the two-point 
function because we assume the non-Gaussianity of the 
CMB fluctuations to be negligible. In addition to space- 
time noncommutativity, the direction-dependent primor- 
dial power spectrum can also be a consequence of many 
phenomena in the early universe, for example see [l9[. 
To compare theoretical prediction with the observation, 
it is convenient to expand the temperature anisotropy 
into spherical harmonics 



AT{n) = T{h) - To = y^^aimYim{n) 



(18) 



where To is the mean temperature of the CMB. If we 
write the direction-dependent primordial power spectrum 
as fil 



r{k) = A{k) 



l + Y.VLMYLM{k) 



LM 



(19) 



the covariancc matrix of aim will be 

LM 

Here, C;^ is the CMB angular power spectrum, given by 

/"OO 

= (47r)M dkeA{k)[Ti,{k)]\ (21) 
Jo 

where Ti-^ is the CMB transfer function which is taken 
to be rotationally invariant. It is known that if the two- 
point function of the temperature anisotropy is not rota- 
tionally invariant, the covariancc matrix of a;,„ will not 
be diagonal. The off diagonal elements of the covariancc 
matrix appear in the second term of eq. (j20p This term 
is given by 



A";"' = (47r)2(-z)'i-'^ / dkeA{k)VLMTi,{k)TiAk), 
Jo 

(22) 



rhm.3 



where £, 



hm.3 



are the Clebsch-Gordan co- 



Iimil2m2' ^Iimil2m2 

efficients. Comparing the power spectrum in eq. (|19|) 

with the one in eq. p7p . we find that the non zero 
components of Vlm are 



3 2 



and V2o^^J-^^n', (23) 



We see that the contributions from spacetime noncom- 
mutativity also influence the diagonal elements of the 



covariancc matrix, i.e. they modify the amplitude of Ci. 
Substituting Vlm from eq. (^5]) into eq. ([^ . we obtain 
the covariancc matrix 



™l^i2m2/ 



1 + 



2 \{2h + 3){2h 



l\ — mf 



1 



- 1 



(24) 



2(2/i-l)V (2^1 -3)(2;i + 1) 



+ ^1 



where Ci,i, = (inf dkk^A{k)Ti,{k)Ti,{k)/Ci,. 

We next constrain the contributions from spacetime 
noncommutativity in this covariancc matrix using the 
five-year WMAP foreground-reduced maps [l^]. Since 
the spacetime noncommutativity also influences the di- 
agonal elements of the covariancc matrix and its con- 
tribution depends on H'^, we write the noncommutative 
contribution in terms of the amplitude of the primordial 
power spectrum ^4^ as 



(25) 



where E = Sn'^rripe'^clL'^/ld. We adopt the procedures 
in IJTj] to compute the posterior probability of the pa- 
rameters E and A^, given the observed temperature 
anisotropics a. 



P(CT|a) cx L{ai\a)P{a), 



(26) 



where a — {Yi,As} is the set of parameters, L{a.\a) is 
the likelihood and P{(j) is the prior. Since the galactic 
contamination cannot be completely removed from some 
regions of the sky, one does not have the full-sky CMB 
maps with well-defined error properties. To reduce the 
galactic contamination, one masks the contaminated re- 
gions as Ci = Mi/STi, where i is the pixel index, Ci is the 
masked CMB map, AT^ is the full-sky map and Mi is 
a mask which is zero at the contaminated points and is 
one elsewhere. This relation can be written in harmonic 
space as 



where the matrix Mi, 



Mir. 



(27) 



LM 



MLM^tkl'm' 



Here, 



/ is given by Mim,i'm' = 
Mim and b im are the 



spherical harmonic coefficients of q. Mi and ATi respec- 
tively. Moreover, due to the instrument noise, the finite 
beam resolution and the discreteness of the temperature 
maps, the contributions to the unmasked CMB map come 
from the sum of the instrument noise with the convolu- 
tion between the signal of the CMB anisotropies and the 
window function, such that 



AT ^Wa + U, 



(28) 



where W is the window function and is the instru- 
ment noise. We suppose that the contributions from the 
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beam and the pixel asymmetries are neghgible. Using 
eqs. (P7|) and the covariance matrix of the masked 
temperature muhipoles can be written as 



where Clm,l'M' = (^^la/^^l'M')- Therefore, the pos- 
terior probabihty function takes the form 



hmil2m2 

+Nlj^mid2 



Pia\A 



LM) OC 



exp I 



det 1/2C 



(35) 



M, 



127712 jl'm' 



(29) 



where Ni^rni,i2m2 is the pixel noise covariance matrix, 
given by iVi^^^,;,™, = AaJ2LM-Lm2m2^LM- Here, 
Aa is the area of each pixel in the temperature map, 
and Nlm is defined as Nlm = J2i ^a{(jl/nf"')YLM{fi), 
where cto is the rms noise of a single observation and nf^^ 
is the number of observations of pixel i. 

In order to compute the likelihood, the inversion of the 
covariance matrix is required. Since the inversion of the 
large matrix is time consuming, we avoid to inverse the 
large covariance matrix by writing the likelihood function 
in terms of the reduced bipolar coefficients instead of Qm- 
The reduced bipolar coefficients are defined as 



Alm — 



E 

Itmtl2m2 



1, "hnu"i2m2^hmii 



(30) 



where dim are the harmonic coefficients of the tempera- 
ture anisotropics. For the full-sky and noiseless case, the 
mean of the reduced bipolar coefficients for the noncom- 
mutative k-inflation can be computed using eq. (j24[) , and 
the non zero components are 



(^oo> 
(A20) 



(31) 



where = Elml'M' ^^lm {C ^) lm,l'M' ""^l'w 

Here, we have used a flat prior on E and a Gaussian prior 
on As- For the Gaussian prior, the mean and variance 
of As are taken from the five-year WMAP results [HI. 
Since the contributions from spacetime noncommutativ- 
ity mainly appear in the (Aqo) and (A20), we restrict the 
multipole index L of to be less than 4. We note that 
the parallel computing of covariance matrices Ci-^rni,i2m2 
and Clm l'M' can be easily implemented in the Healpix 
package 0|. The data that are used to computed 5 Alm 
are the five-year WMAP foreground-reduced V2 and Wl 
differential assembly temperature maps (20| . These maps 
are masked using the band-limited masks in [ll[ . Accord- 
ing to we limit the multipole index of c;™ and M^m 
to be ^ < 62 and L < 92 respectively. The covariance ma- 
trix C can be computed from the covariance matrix in eq. 
(|24p . Since the spacetime noncommutativity does not af- 
fect the cosmic evolution after the inflationary epoch, we 
can use CMBEASY [23| to compute the CMB transfer 
function and C; by supposing that the cosmic evolution 
after inflation obeys the ACDM model whose parameters 
are taken from the best fit value of the five-year WMAP 
results [2l|. However, recall that the amplitude of the 
primordial power spectrum is treated as a free parame- 
ter. We compute the posterior probability function for 
a and marginalize it over A^ to obtain the marginalized 
posterior probability function for E. The marginalized 
posterior probability functions for E obtained from V2 
and Wl maps have a peak at negative E, which can oc- 



2Al'E,''^{—iyCiA — —. (32) cur if is negative. However, we cannot use the above 

;>4 V -i^(^^ ~ -i) analysis to constrain parameters Ls, and e simultane- 



In the more realistic case, the mask and the instrument 
noise must be taken into account, so that the reduced 
bipolar coefficients are computed from ci„i in eq. ([27]) 
and their mean is 

(Alm) = ^ {~'^)"''^'^hini,l2m2^li^mil2-m2- (^3) 
Iil2'mim2 

It can be seen from the above equation that due to the 
effect of the mask, {Aqq) and {A20) will not be the only 
non zero components of {Alm)- 

We define 5 Alm — Alm — (^la/), and write the co- 
variance matrix for 5 Alm as 



ClKLL'M' — ^ ^ ^hmil2m2^hi 

Iimil2m2 ^37713/47714 



M' 

7773/47774 



Cil777i,;3 7773Ci27772,U7774 + i 777i ,(4 7774 ^'(2 7772 ,i 3 77l3 ) 1 (34) 



analysis to constrain parameters Ls, cj. and e simultane- 
ously due to degeneracy among these parameters. Since 
we are interested in the constraint on Ls, we suppose 
that values of and e are known. For canonical scalar 
field, E > because — 1. Hence, we restrict ourselves 
to the case where E > 0, i.e. > 0. In this case, the 
confidence intervals for E can be obtained from the areas 
under the curves of The marginalized posterior probabil- 
ity functions in figure 1. 

In order to estimate the upper bound for Ls from the 
upper bound for E in table 1, we further assume that 
the values of the product ecs are known precisely, so that 
the upper bound for Ls can be written as < 1.4 x 
104[e2c2m4]-i/4Gev-i ~ 1.4x 10-28[ecs]-i/2cm at 99.7% 
confidence level. As shown in [l^, the ratio of the tensor 
to scalar perturbations amplitudes r depends on eCs . We 
roughly estimate the values of eCs using the values of 
r from the five-year WMAP results, and obtain Lg < 
10^'^^cm at 99.7% confidence level. 
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68%CL 


95%CL 


99.9%CL 


Wl 


0.024 


0.048 


0.070 


V2 


0.029 


0.055 


0.078 



TABLE I: The upper bound for the parameter E x 10 
from Wl and V2 maps. 

7 I 1 1 1 1 1 1 




0.12 0.14 



S X 10-'" 



FIG. 1: The marginahzed posterior probabihty functions for 
E. The solid line represents the probability function from Wl 
map, while the long dashed line represents the probability 
function from V2 map. Here, E is restricted to be positive. 

IV. CONCLUSIONS 

In this work, we study the effect of spacetime noncom- 
mutativity on the rotational invariance of the primordial 
power spectrum, and constrain the contributions from 
this effect using CMB data. In the slow-roU approxi- 
mation, the deviation from rotational invariance of the 



primordial power spectrum due to the spacetime non- 
commutativity effect depends on the factor 3H'^L'^/16, 
but not on sound speed. Here, H is evaluated at the 
time when perturbation mode k crosses the horizon dur- 
ing inflation. 

Since the primordial power spectrum is direction- 
dependent, the covariance matrix for the harmonic coeffi- 
cients of the CMB temperature anisotropics has off diag- 
onal elements. As is well known, this implies that statis- 
tics of the CMB anisotropics become anisotropic. The 
spacetime noncommutativity also contributes to the di- 
agonal elements of the covariant matrix suggestting that 
the noncommutative contribution also modifies the am- 
plitude of the CMB angular power spectrum. 

Both contributions from spacetime noncommutativity 
are simultaneously constrained using five-year WMAP 
foreground-reduced V2 and Wl maps. The upper 
bound for the quantity Ls[ecs]^/^ is approximately 
1.4 X 10~^^cm at 99.7% confidence level. Estimating the 
values of ec^ using the tensor to scalar perturbations 
amplitudes from the five-year WMAP results, we obtain 
Ls < 10"^'' cm at 99.7% confidence level. 
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